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A Method  for  the  Calculation  of 
Parachute  Opening  Force*  for 
High  Altitude  Balloon  Payloads 


I.  INTRODUCTION 


Each  year  a large  number  of  valuable  scientific  payloads  are  committed  to 
flight  on  high  altitude  balloons.  In  general  the  balloons  ascend  into  the  stratosphere. 
They  have  gone  as  high  as  170,  000  ft  MSL  (mean  sea  level)  and  often  have  flight 
durations  of  several  days.  At  the  conclusion  of  the  flight  the  payloads  are  usually 
recovered  by  means  of  a parachute  system  for  later  reuse.  The  parachute  apex 
is  attached  to  the  bottom  end  fitting  of  the  balloon  and  flown  unpacked  In  line  with 
the  payload  suspended  below.  At  the  conclusion  of  the  balloon  flight  the  parachute 
Is  separated  from  the  balloon  and  descends  with  the  payload  to  the  ground. 

One  of  the  questions  most  frequently  asked  by  the  sclentists/experlmenters 
who  design  the  payloads  and  Interpret  the  data  obtained  concerns  the  magnitude  of 
the  parachute  opening  force.  Conventional  methods  for  deriving  parachute  per- 
formance parameters  in  the  dense  lower  atmosphere  do  not  reliably  predict  this 
information  for  the  parachute  released  in  the  very  high  altitude  environment. 

An  analytical  answer  to  the  question  is  sought.  Existing  theory  and  an  analytical 
method  recently  developed  by  Ludtke1  for  horizontal  deployment  is  adapted  to  the 
balloon  system  requirements.  Engineering  methods  are  used  to  develop  expres- 
sions which  allow  the  assessment  of  the  parachute  opening  forces. 

(Received  for  publication  7 May  1975) 

1.  Ludtke,  W.P.  (19 1 2)  A Technique  for  the  Calculation  of  the  Opening-Shock 
Forces  for  Severs f *** r s n-n  « i ac 


2.  OPENINU  DYN\MIC:s 


At  the  termination  of  the  flight  when  the  apex  o.  the  parachute  system  is 
separated  from  the  balloon,  the  parachute/payload  system  begins  to  free  fall.  The 
strain  energy,  which  had  been  stored  in  the  parachute  suspension  lines,  riders,  and 
canopy  when  the  system  was  supported  by  the  balloon,  causes  the  parachute  to  be 
accelerated  towards  the  payload.  As  the  system  descent  velocity  increases  a few 
parachute  gores  which  ma>  begin  to  catch  air  and/or  skin  drag  on  the  canopy  tend 
to  stretch  the  system  back  into  line.  The  payload  will  see  a momentary  force  peak 
as  the  system  snaps  back  into  tension.  This  "snatch"  force  is  generally  taken  to  be 
the  beginning  of  parachute  filling. 

The  further  motion  of  the  descending  system  can  then  be  described  by  Newton's 
second  law  of  motion. 


EF  * ma  , 


l 2 
- JO 


W dv 
g dt  * 


(1) 


Integration  of  Eq.  (1)  requires  a knowledge  of  the  change  of  air  density  (p)  with 
time  and  of  the  change  of  drag  area  (CpS)  with  time.  Let  it  be  assumed  that  during 
the  opening  period  the  atmospheric  density  can  be  represented  by 


where  a j and  a2  are  constants  and  z is  altitude  MSL.  Then  altitude  at  any  instant 
may  be  determined  from 


dz 

dt 


v(t)  . 


(3) 


3.  LI' DIKE'S  TECHNIQUE 

3. 1 Opening  Force  Ratio 

Ludtke1  has  suggested  that  parachute  opening  can  be  considered  with  respect 
to  a reference  opening  time,  tQ.  This  reference  time  is  the  instant  when  the 
canopy  is  considered  to  have  just  achieved  its  fully  inflated  steady-state  diameter 
and  volume  for  the  first  time.  Additional  inflation  beyond  the  reference  time  is 
considered  to  be  caused  by  the  elongation  of  the  canopy  materials  under  the  applied 
loads.  For  purposes  of  curve  fitting  the  reference  time  (tQ)  was  established  as 
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the  first  point  on  the  infinite  mass  wind  tunnel  force-time  plot  where  the  insta  >tan- 
eous  force  (F)  was  equal  to  the  steady-state  drag  force  (Fs). 


F =ivv2rDs  , 


F.^pv2  CdSq  . 


i 

A 


For  infinite  mass  conditions  where  p is  constant  and  v = vg  it  was  determine  ! that 


F _ 


cL  S 
Vo 


(4) 


for  most  solid  parachutes. 

Furthermore,  all  parachutes  have  some  initial  drag  area  at  line  stretch 
(v  : v3).  Therefore  it  is  suggested  that 


R <n>  * ^(l-n)^f-^  + nj 


> (5) 


where  p is  defined  as  the  ratio  of  the  projected  mouth  diameter  (D^)  at  line  st  etch 
to  the  steady-state  frontal  diameter  (DMq).  Expanding  Eq.  (5) 


R(n)  = (1-n) 


+ 


2rj  (1-rj) 


+ n 


(C) 


This  expression  has  been  shown  to  be  valid  for  the  finite  mass  case  as  well 


1.2 


3.2  Determination  of  Reference  Time 

At  any  given  time  during  the  parachute  uflation  the  parachute  drag  force  is 
proportional  to  the  square  of  the  maximum  inflated  diameter.  Using  this  observa 
tion  the  following  assumptions  can  be  made, 

(a)  The  ratio  of  the  instantaneous  mouth  inlet  area  (A^j)  to  the  steady-state, 
fully  Inflated  mouth  area  (Aj^0)  is  in  the  same  proportion  as  the  instantaneous  dra  j 
area  ratio. 


r^-  = R(r>)  . (7) 

Mo 


i 


2. 


Berndt,  R.J.,  and  DeWesse,  J.H.  (1966)  Filling  time  prediction  app; •>  ch 
for  solid  cloth  type  parachute  canopies,  AIAA  Aerodynamic  Deceleration 
Systems  Conference,  Houston,  Texas. 
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(b)  The  ratio  of  the  instantaneous  pressurised  cloth  surface  (Ag)  to  the  canopy 
surface  area  (AgQ)  is  in  the  same  proportion  as  the  instantaneous  drag  area  ratio. 

r^-  * R(n)  . (8) 

*So 

(c)  The  rest  of  the  canopy  Is  assumed  to  be  unpreasurized  and  therefore  has 
no  net  air  flow. 

Applying  the  principle  of  conservation  of  mass 
dm  - - mou|  . 

pa?  =pvAMoR(n) -pAsoR(',,k(¥)  y2n  • <9> 

where  V is  the  instantaneous  volume  of  air  captured  by  the  parachute  and  k,  n,  and 
Cpare  coefficients  which  are  a function  of  the  canopy  cloth  permeability.  A more 
detailed  explanation  of  these  coefficients  may  be  found  in  reference  1. 


4.  CALCULATION  OF  REFERENCE  TIME 


Parachute  opening  data  suggests  that  the  initial  effective  diameter  ratios  used 
in  Eqs,  (7)  and  (8)  may  not  be  the  same  as  the  initial  drag  area  ratio.  Equations 
(7)  and  (8)  then  become 


R(V 


(10) 


and 


•Ac 

j—  » R(hs)  , (11) 


Equation  (9)  then  becomes 


P dlf  P v^Mo  R ~ p ^So  R^S^  y2n  ‘ 


8 


The  reference  time  may  now  be  calculated  by  making  an  initial  guess  for  t 
and  numerically  integrating  Eqa.  < 1),  (3),  and  (12)  simultaneously  from  t - 0 
and  V * 0 to  V •*  VQ.  The  quantity  is  the  known  steady  ^tate  parachute  volume 
which  may  te  determined  . shown  in  reference  3.  Then  the  time  V * VQ  la  a new 
estimate  for  and  the  equations  can  again  be  integrated  to  obtain  a new  estimate 
for  tQ.  This  process  will  converge  to  tQ. 


3.  END  POINT  DETERMINATION 


Once  the  parachute  reaches  its  nominal  shape  it  continues  to  increase  in  area 
according  to  Eq.  (4).  This  inflation  process  is  limited  by  the  loads  applied  to  the 
canopy,  the  elasticity  of  the  canopy,  and  its  constructed  strength  (Fc).  A linear 
load  elongation  relationship  is  used  to  determine  the  maximum  drag  area  at  any 
time. 


E emax 

F ' F 

c 


(13) 


where  e is  the  strain  of  the  parachute  material  when  loaded  to  its  constructed 

max  r 

strength  (Fc).  At  any  time  the  force  (F)  on  the  parachute  is  given  by  the  drag. 

F=|pR(n)CDSov2  , (14) 


In  steady  descent  the  applied  load  is  just  the  system  weight  CW)  and  the  para- 
chute drag  area  is  increased  by  the  factor 


V! 


cmax 

F 

c 


However,  it  must  be  recalled  tb3t  parachute  drag  coefficients  are  usually  deter- 
mined using  the  parachute  nominal  drag  area.  Therefore,  at  any  given  time  the 
drag  area  ratio  is  limited  to 


3.  Ludke,  W,  P.  (1970)  A New  Approach  to  the  Determination  of  the  Steady -State 
Inflated  Shape  and  Included  volume  of  Several  Parachute  Types  in  24-Gore 
~ 6ore  Configurations,  NOLTR  70-178, 


9 


t 


c 


The  drag  area  ratio  of  a parachute  thus  initially  increases  according  to 
Kq.  (6)  until  t = t and  then  according  to  Eq.  (4).  The  increase  in  drag  area  ratio 
is  limited  to  the  value  given  by  Eq,  (15)  at  any  instant  of  time.  The  calculation  of 
the  limiting  drag  area  ratio,  therefore,  requires  the  simultaneous  solution  of 
Eqs.  (14)  and  (15). 

6.  \PPI.IC4TION8 

In  order  to  evaluate  the  behavior  of  these  equations  a computer  program  was 
written  to  numerically  integrate  the  equations  as  a function  of  time.  Table  1 shows 
the  input  parameters  required  for  a 135-ft  diameter  flat  circular  parachute. 

Figures  1 through  5 show  the  force,  velocity,  and  drag  area  ratio  as  a function  of 
time  for  the  135-ft  parachute  opening  at  a snatch  velocity -<v  ) of  160  ft/sec  in 
vertical  deployment  at  several  snatch  altitudes.  It  can  be  observed  that  the  system 
velocity  increases  after  snatch  until  the  drag  force  exceeds  tne  system  weight. 

The  maximum  force  (F„  ) occurs  well  before  the  parachute  is  fully  open  at  low 

altitudes.  At  1000  ft  M3L,  Fmax  occurs  when  the  drag  area  is  only  10  percent  of 


Table  1.  Input  Parameters  for  135-ft  Parachute 


Do 

135  ft 

W 

10,  000  lbf 

CD 

1.22 

vo 

290,  000  ft" 

F 

88.  000  lbf 

c 

£ 

0.30 

max 

"s 

0.02 

'’m 

0.02 

0.674 

CP 

1.15 

k 

1.46042 

n 

0.63246 
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f 
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Figure  5.  Opening  of  a 135-ft  Diameter  Parachute  at  100000  ft  MSL 


of  the  nominal  drag  area.  As  the  snatch  altitude  is  inc-eased,  however,  the  time 

of  maximum  force  occurrence  becomes  later  and  is  after  the  reference  time  (t  ) 

o 

at  the  higher  altitudes. 

Figure  6 illustrates  that  for  the  135-ft  parachute  the  reference  filling  time  (tQ) 
decreases  with  increasing  altitude.  Above  approximately  60,  000  ft  MSL  the  filling 
time  is  essentially  constant. 

In  Figure  7 the  maximum  forces  calculated  during  opening  of  the  135-ft  para- 
chute in  the  vertical  and  horizontal  modes  are  compared.  The  difference  between 
the  modes  of  deployme...  nay  be  partially  explained  by  the  fact  that  at  the  end  of 
the  opening  calculation  (t  » tQ)  the  drag  fcrce  tends  towards  W for  the  vertical 
deployment.  The  maximum  opening  force  as  a function  of  altitude  has  a peak  at 
approximately  55,  000  ft.  The  maximum  force  is  limited  by  the  maximum  drag 
area  and  velocity  experienced.  In  horizontal  deployment  the  maximum  drag  force 
is  approximately  F , which  is  the  force  that  would  be  experienced  by  the  fully 
open  parachute  at  v . For  vertical  deployment  F exceeds  F , due  to  the  in- 
crease  in  velocity  after  snatch.  The  maximum  force  in  horizontal  or  vertical  de- 
ployment may  also  exceed  Fg  due  to  stretching  of  the  parachute  (R  > 1), 

As  a second  example  consider  the  input  parameters  shown  in  Table  2 for  a 
35-ft  parachute.  Figures  8 and  9 show  Fmax  for  the  opening  of  the  3 5 -ft  parachute 
in  horizontal  deployment  assuming  a constant  initial  velocity  (v  r 400  ft /sec) 
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Figure  6.  Release  Time  (tD) 
as  a Function  of  Altitude  for 
a 13  5-ft  Diameter  Parachute 


Figure  7.  Maximum 
Opening  Force  for  a 3 5-ft 
Diameter  Parachute  in 
Horizontal  Deployment 
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Figure  9.  Maximum  Open  Force  for  a 35-ft  Diameter  Parachute 
in  Horizontal  Deployment 


o 

and  a constant  initial  dynamic  pressure  (q  = 190  lbf/ft  ).  Note  that  Fmftx 
exceeds  F for  the  v = 400  ft/sec  case  due  to  the  straining  induced  in  the 

3 8 o | 

parachute  canopy.  Also  r.^te  that  the  maximum  force  for  the  q = 190  lbf/ft* 

example  continues  to  iticrease  with  altitude  and  is  in  agreement  with  the 

4 

trends  discussed  elsewhere  in  the  literature. 

A number  of  cases  were  evaluated  at  different  initial  velocities  and  alti- 
tudes for  the  135-ft  parachute  and  are  plotted  in  Figure  10.  The  maximum 
spread  in  the  Fmax  data  occurs  in  the  region  of  maximum  value.  If  the  data 

of  Figure  10  are  replotted,  using  a free  fall  time  (t  ) to  calculate  v , then 

s s 

the  result  is  as  shown  in  Figure  11. 


4.  Performance  nf  Design  Criteria  for  Deployable  Aerodynamic  Decelerators. 
(1963)  TP  ASi:-TR-6i-579,  AFFDL,  AFSC. 
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7.  CONCLUSIONS 


If  U can  be  assumed  that  the  snatch  force  will  occur  at  a relatively  constant 
time  alter  separation  from  the  balloon  then  the  plot  of  Figure  11  shows  that  it  would 
be  better  to  terminate  the  mission  at  higher  altitudes  rather  than  at  lower  altitudes. 
Although  it  is  not  known  whether  the  snatch  time  (tg)  is  constant  with  altitude.  It 
is  known,  in  general,  to  be  less  than  5 sec  based  on  limited  flight  data  et.  altitudes 
above  approximately  50,  000  ft  MSL. 

No  attempt  is  made  to  quantitatively  prove  or  disprove  the  validity  of  the 
mathematical  model  herein  discussed.  It  does  qualitatively  agree  with  experience 
and  its  quantitative  results  are  a function  of  the  input  parameters,  which  may  not 
all  be  well  known.  However,  the  trend  analysis  provided  by  the  type  of  plot  shewn 
in  Figures  10  and  11  is  useful  for  high  altitude  balloon  mission  planning. 
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Symbols 

o 

Instantaneous  pressurized  canopy  surface  area,  ft 

2 

Steady-state  inflaved  canopy  surface  area,  ft 

2 

Instantaneous  canopy  mouth  area,  ft 

2 

Steady-state  inflated  mouth  area,  ft 
2 

Acceleration,  ft/sec 
Constants,  Eq.  (2) 

Parachute  coefficient  of  drag 

Parachute  pressure  coefficient;  relates  internal  and  external  pressure 
(P)  on  canopy  surface  to  the  dynamic  pressure  of  the  free  stream 

2 

Instantaneous  drag  area,  ft 

2 

Steady-state  drag  area,  ft 

Instantaneous  diameter  of  the  mouth  of  the  canopy,  ft 

Steady-state  inflated  mouth  diameter,  ft 

Nominal  parachute  diameter,  ft 

Instantaneous  force,  lbf 

Constructed  strength  of  the  parachute,  lbf 

Maximum  opening-shock  force,  lbf 
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Steady-state  drag  force  that  would  be  produced  by  the  fully  open 
parachute  at  v , Ibf 

2 

Gravitational  acceleration,  ft/sec 
Permeability  constant  of  canopy  cloth 
Mass,  slugs 

Permeability  constant  of  canopy  cloth 

o 

Dynamic  pressure,  Ibf/ft 

Ratio  of  instantaneous  drag  area  to  nominal  drag  area  (Cj-^) 

2 

Instantaneous  inflated  canopy  surface  area,  ft 

2 

Canopy  surface  area,  ft 
Instantaneous  time,  sec 

Reference  time  when  the  inflating  parachute  has  reached  the  design 
drag  area  for  the  first  time,  sec 

Fall  time  before  snatch,  sec 

3 

Instantaneous  volume  of  air  collected  by  canopy,  ft 

Total  volume  of  air  which  must  be  collected  during  the  inflation 
process,  ft3 

Instantaneous  system  velocity,  ft/sec 
Snatch  velocity,  ft/sec 
System  weight  j( 

Altitude,  ft  M3L 


Instantaneous  strain 

Maximum  strain 

Ratio  of  parachute  projected  mouth  diameter  at  line  stretch  to  the 
steady-state  inflated  mouth  diameter 

Ratio  of  initial  mouth  diameter  to  steady-state  inflated  mouth 
diameter 

Ratio  of  initial  effective  pressurized  cloth  diameter  to  canopy  nominal 
diameter 

3 

Air  density,  slugs /ft 
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